An Improved Standard Model Prediction Of BR(B tv) 
And Its Implications For New Physics 
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The recently measured B ^ tv branching ratio allows to test the Standard Model by probing 
virtual effects of new heavy particles, such as a charged Higgs boson. The accuracy of the test 
is currently limited by the experimental error on BR{B — >■ tv) and by the uncertainty on the 
parameters fs and |Kib|- The redundancy of the Unitarity Triangle fit allows to reduce the error 
on these parameters and thus to perform a more precise test of the Standard Model. Using the 
current experimental inputs, we obtain BR{B — >■ tv)sm = (0.84 ± 0.11) x 10"**, to be compared 
with BR{B -5> rj/)exp = (1.73 ± 0.34) x 10~*. The Standard Model prediction can be modified by 
New Physics effects in the decay amplitude as well as in the Unitarity Triangle fit. We discuss how 
to disentangle the two possible contributions in the case of minimal fiavour violation at large tan /3 
and generic loop-mediated New Physics. We also consider two specific models with minimal fiavour 
violation: the Type-II Two Higgs Doublet Model and the Minimal Supersymmetric Standard Model. 



INTRODUCTION 



Flavour physics offers the opportunity to probe vir- 
tual effects of new heavy particles using low-energy phe- 
nomena, involving Standard Model (SM) particles as ex- 
ternal states. New Physics (NP) can generate large ef- 
fects in Flavour Changing Neutral Currents (FCNC) and 
CP violating phenomena even for NP particle masses 
much above the TeV scale, if new sources of flavour and 
CP violation besides the Yukawa couplings are present. 
The strong NP sensitivity is mainly due to the Glashow- 
Iliopoulos-Maiani (GIM) suppression of FCNC processes 
in the SM [T]. However, other suppression mechanisms 
can be at work in the SM, making a few non-FCNC de- 
cays interesting for NP searches. In particular, the he- 
licity suppression of the charged current decay B ^ tv 
makes it potentially sensitive to the tree-level effects of 
new scalar particles |2j . A typical example is given by the 
exchange of charged Higgs bosons in multi-Higgs exten- 
sions of the SM, such as the type-H Two Higgs Doublet 
Model (2HDM-n) or the Minimal Supersymmetric Stan- 
dard Model (MSSM), in the large tan/3 regime. 

In the SM, the branching ratio oi B tv can be 



written as: 

BR{B ^ TV) = [''^) ^^'^""1 • 

(1) 

The Fermi constant Gp, the B (r) mass tob i'm-r) and 
the B lifetime tb are precisely measured [3]. The de- 
cay constant of the B meson /b is known with C(10%) 
uncertainty. We use the lattice QCD (LQCD) average 
Jb = 200 ± 20 MeV g]. Concerning the error attached 
to lattice averages, we combine in quadrature the statis- 
tical and systematic errors, assuming Gaussian distribu- 
tions. This is justified since present lattice systematic 
errors arise from the combination of several independent 
sources of uncertainty. Therefore they are well described 
by a Gaussian distribution, no matter what the distribu- 
tions of the individual sources are. ^ 

The absolute value of the Cabibbo-Kobayashi- 
Maskawa (CKM) [5] matrix element Vub is determined 
from the measurements of the branching ratios of ei- 



Notice that in the past we used to assign a flat distribution to 
the lattice systematic errors, since they were dominated by the 
uncertainty associated to the quenched approximation. 
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elusive and inclusive semileptonic b ^ u decays. Its 
precision is limited by the uncertainty of the theoreti- 
cal calculations. Although inclusive determinations are 
systematically higher than exclusive ones, the two values 
are compatible, once the spread of inclusive determina- 
tions using different theoretical models is considered. For 
the exclusive decays, we use the HFAG averages E] 

BR{B 7r£i/),2<i6GcV2 = (0-94 ± 0.05 ± 0.04) x 10"'' , 



BR{B tt£v] 



g2>l6GcV2 



(0.37 ±0.03 ±0.02) X 10" 



together with the theoretical estimates of the relevant 
normalized form factors 

FF(g2 < leGeV^) = 5.44 ± 1.43 [S] , 
FF(g2 > 16 GeV^) 2.04 ± 0.40 g] , 

to obtain iKbT'"''' = (33.3±2.7) x IQ-^. For inclusive de- 
cays, we quote |Kb|'"''' = (40.0± 1.5±4.0) x 10"", where 
we define the second error as a flat range accounting for 
the spread of the different models [9j. 

Our grand average of inclusive and exclusive determi- 
nations is IK&I = (36.7 ± 2.1) x 10-*, obtained fr om the 
probability density function (p.d.f.) in Fig. [l] From this 
p.d.f. we get 



BR{B Tv) = (0.98 ± 0.24) x 10" 



(2) 
tlO) at 



compatible with Bi?exp = (1-73 ± 0.34) x 10" 
~ I.Bct. 

A few percent precision is expected to be reached by 
LQCD using Petaflop CPUs for fs and the form factors 
entering the exclusive determination of \Vub\ HJ- Con- 
sidering how challenging the measurement of BR(B — > 
Tif) in a hadronic environment is, it is difficult to imag- 
ine a similar improvement in precision of the experimen- 
tal measurement, unless a SuperB factory will be built, 
leading also to a better direct determination of \Vub\ [H]. 
On the other hand, it has been pointed out in Ref. [12] 
that the indirect determination of \Vub\ from the Unitar- 
ity Triangle (UT) fit in the SM is more accurate than the 
measurements, yielding a central value close to the exclu- 
sive determination. Therefore a more precise prediction 
of BR{B Tv) in the SM can be obtained combining 
the direct knowledge of \Vub\ and /s with the indirect 
determination from the rest of the UT fit. 



UTFIT-IMPROVED STANDARD MODEL 
PREDICTION 

In the UT fit [HITl], CP-conserving and CP-violating 
measurements are combined to constrain p and f]. The fit 
also provides an a-posteriori determination of \Vub\ which 
includes the direct measurement as well as the indirect 
determination from the other constraints. Similarly, an 
improved determination of fs from both LQCD and ex- 
perimental constraints is obtained \X2i- 



= 0.002 



0) 



>■ 0.001 5 



(0 



0.001 



0.0005 




0.003 0.0035 0.004 0.0045 0.005 



|V.J 



ub' 



FIG. 1: P.d.f. of \Vub\ obtained combining inclusive and ex- 
clusive measurements of the b ^ u semileptonic decays. The 
dark (light) region corresponds to the 68% (95%) probability 
interval. 



The most accurate prediction of BR{B — > riy) in 
the SM can then be obtained performing the SM fit 
without including the measurement of BR{B tv) 
as a constraint. The fit gives p ~ 0.149 ±0.021 and 
f] = 0.334 ± 0.013 together with fs = (196 ± 11) MeV 
and iKbl = (35.2±1.1) x lO""*. The posterior p.d.f.'s are 
shown in Fig. [2] 

The same SM fit gives the p.d.f. in Fig. |3j from which 
we obtain 



BR{B Tv)sM = (0.84 ± 0.11) x 10"^ 



(3) 



In Fig. |4]we present the compatibility plot for BR{B — > 
'''^)sM- The colored regions represent the pull from the 
UT fit result. The present experimental value, repre- 
sented by a cross in the plot, displays a deviation of 
~ 2.5(7. This deviation can be interpreted as a similar 
same-sign statistical fluctuation (or a correlated system- 
atic error) in BaBar and Belle results or as a hint of NP 
effects. A more definite answer needs new data to be 
collected. 

From Eq. ([s]), one can easily predict the SM value of 
BR{B p,y) and BR{B ev). We obtain 



BR{B 
BR{B 



■ ^iy)su = (3.8 ±0.5) X 10" 
. ev)sM = (8.8 ± 1.2) X 10 



(4) 
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The precision on the experimental measurements |15j is 
still far from probing such small values. The current best 
limits are BR(B fiv) < 1.0 x 10"^ [T^ and BRiB 
eiy) < 1.0 X 10-6 [7] at 90% C.L. 
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FIG. 2: Posterior p.d.f. for \Vub\ (top) and /s (bottom), ob- 
tained from the UT fit, without taking BR{B — >■ tv) as input. 
The dark (light) region corresponds to the 68% (95%) proba- 
bility interval. 



FIG. 3: P.d.f. for BR{B -s- riy) predicted using the UT fit. 
The dark (light) region corresponds to the 68% (95%) proba- 
bility interval. 




BR(B^xv)[10-1 



MODEL-INDEPENDENT PREDICTIONS 

Let us assume in the following that NP is at work. In 
this case, the prediction in Eq. ^ could be modified by 
i) NP effects in the decay amplitude and/or ii) NP ef- 
fects in the UT fit. If more precise measurements will 
provide evidence of a discrepancy, one should be careful 
in interpreting it as evidence of NP in the B tv de- 
cay amplitude. In fact, other inputs of the UT analysis 
(for example Am^ {q — d,s)) might be affected by the 
presence of contributions beyond the SM. We would like 
to disentangle the two possible NP effects. To this aim, 
we compute the prediction of BR{B — > tv) in several NP 
scenarios assuming that NP contributions to the B ^ tv 
decay amplitude are negligible. This prediction will be 
denoted as i?i?modoi- A discrepancy between BR^^dci 
and BRexp would unambiguously reveal NP contribu- 
tions to the B ^ TV decay amplitude in the considered 



FIG. 4: Compatibility plot for BR{B tv). The cross 
marks the current world average; colours give the agreement 
(in number of a) with the data-driven SM prediction. 



scenario. 

As is common practice in the literature, we also provide 
results in terms of the ratio 



pcxp BRcxp 

^modcl 



BR 



(5) 



model 



The use of -Rmodci particularly convenient for NP mod- 
els with Minimal Flavour Violation (MFV) [I71[T5], de- 
fined as models where the only source of flavour violation 
are the quark masses and the CKM matrix [TH]. Indeed, 
BRmfVj the full prediction of the branching ratio in- 
cluding NP in the decay amplitude, and BRmfv have 
the same dependence on \Vub\ and Jb, so that they can- 
cel in the ratio i?MFV = BRmfy/BRmfy- Therefore, 
Rmfv can be computed theoretically without specifying 
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FIG. 5: P.d.f. of R^^rj, obtained using the UUT construction. 



the value of | Vub \ and /b ■ Rmfv is constrained by -Rmfv ' 
which contains the experimental error as well as the un- 
certainty on \Vub\ and /b- 

Following Ref . , we distinguish several scenarios ac- 
cording to the NP flavor structure. In each scenario, we 
remove all the inputs that might be affected by NP from 
the UTfit-based determination of BR{B — > tv). This 
gives a NP-independent prediction of BR. 

In MFV models one expects the tree-level processes 
and the angles of the UT not to deviate from the SM 
prediction, while the values of Am^ and ex are expected 
to change.^ We can then replace the full SM UT fit 
with the Universal UT (UUT) construction [21]. In 
the case of the UUT, the knowledge of /b is given by 
LQCD only, resulting in a larger error on i3i?uuT- Us- 
ing the currently available experimental inputs, we ob- 
tain i?i?uuT = (0.87 ± 0.20) X 10"'' corresponding to 
-^UUT ~ ^ '^•6' shown in Fig. [s] (for comparison, 
see the SM result in Eq. ([s])). Clearly, the determination 
of BRtjut will benefit considerably from the expected 
improvements in future LQCD calculations. 

In MFV models with one Higgs doublet (or two Higgs 
doublets at small tan/3), one expects negligible NP ef- 
fects in the B ^ tu decay amplitude, while a deviation 



could be induced on Am^, Atos, and ck- Should -R™ut 
deviate from one significantly, these models would then 
be excluded. 

In the case of MFV models with two Higgs doublets 
at large tan/3, the value of i?^ux could be shifted from 
one by the contribution of the charged Higgs boson to 
the decay amplitude. 



^ In MFV models one has to assume that the large measured value 
of the Ba mixing phase is a statistical fluctuation. Otherwise, 
MFV would be excluded |20) . 



CONSTRAINTS ON 2HDM-II 

As an explicit example of the discussion above, we con- 
sider the 2HDM-II. In this model, the interaction be- 
tween quarks and the charged Higgs is defined by 
the Lagrangian 



3 

^—^ Vtanri 



/3--"' 2 



-h tan ^ md, ^-^^^ + H.c. 

and FCNC are absent at the tree level. 
We can write 1251: 



(6) 



-R2HDM 



1-tanV- 



(7) 



where is the mass of the charged Higgs boson. 

Eq. ([7]), together with the p.d.f. of -R^ux provided by 
the UUT fit, gives a constraint on tan/3/TO^+ as shown 
in Fig. [6j The charged Higgs contribution typically sup- 
presses BR[B — >■ Tv) with respect to the SM, contrary to 
current experimental results. An excess can be obtained 
if tan /? > \f2 'rnfj+ /m b (corresponding to the rightmost 
peak in the left plot of Fig. [6) tan/3 = (29±2) mjy+/(100 
GeV)), yielding an upper limit on to//+ for a given value 
of tan/3. The current direct searches give a lower 
limit of mH+ > 79 GeV at 95% C.L. [S], while the 
measurement of BR{B — ^ Xg'j) imphes mH+ > 295 
GeV at 95% C.L. for the 2HDM-II charged Higgs bo- 
son [21]. This bound excludes the rightmost peak in 
Fig. [6] for tan /3 < 80. In addition, one can consider the 
bound on tfm (3 /mH+ from BR{B ->■ Dtu)/BR{B 
Dtu) where £ denotes light leptons [2 |25]. Using 
the world average (49 ± 10)% [2^ and formula (9) of 
Ref. 23 we obtain the following 95% probability re- 
gions for tan/3/mjj+: tan/3/mjj+ < 0.17 GeV""^ and 
0.46 GeV" 1 < tan/3/mff+ < 0.55 GeV"^ (see the right 
plot in Fig. [6j In this for the B TV bound, 

there is an allowed region at large tan/3/m//+. Assuming 
flat priors in [5,120] for tan/3 [28] and [100,1000] GeV 
for TTijLf + , we obtain the plot in Fig. [t] For tan /3 > 22 
B ^ TV gives a lower bound on tojj+ stronger than the 
one from B — > Xg^- The fine-tuned regions for large 
ta,n/3/mB+ allowed individually by the B tv and the 
B — >■ Dtv constraints do not overlap and are therefore 
excluded. We thus obtain an absolute bound 



tan 13 < 7.4 



100 GeV 



(8) 



In addition, we compute the prediction for BR{Bs 
/i+/x~) and obtain 

BR{B, fi+fi-) = (4.3 ± 0.9) X 10"^ 



(9) 



([2.5,6.2] X 10"^ @95% prob.). 
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FIG. 6: P.d.f. of tan/3/m^+ computed from Eq. Q and the fit result for i^^ (left), or using BR{B Dtv)/BR{B Dlh 
(right). 



Our results are in agreement with Ref. |29l , where 
the effect of BR{B — >■ ri/) and other constraints on the 
2HDM-II has been recently analysed. However, our anal- 
ysis differs in several aspects. First, in the UUT analysis 
we keep all the angles, which are unaffected by MFV 
NP effects.'^ Second, we neglect sub-percent contribu- 
tions to tree-level decays, allowing us to use all determi- 
nations of \Vub\- Third, we only consider the dominant 
constraints from B — >■ and B — >■ tv. Finally, we use 
the Bayesian approach detailed in Ref. [5T| . 

One of the most interesting features of the relation in 
Eq. ([7| is that it does not depend on the flavour of the 
final lepton [22] , since the helicity suppression in the SM 
compensates the scaling of the Higgs couplings with the 
mass. This means that, provided the evidence of a dis- 
crepancy in B — > Tiy, the same effect should be observed 
in i? -> {i — e, ^). For these decays, we get 



scenario 


IKbl X 10* 


fB (MeV) 


BR X 10* 


pull 


UT 


35.2 ± 1.1 


196 ± 11 


0.84 ± 0.11 


2.5o- 


UUT 


35.0 ± 1.2 


200 ± 20 


0.87 ±0.20 


2.2cr 


no-fit 


36.7± 2.1 


200 ± 20 


0.98 ± 0.24 


1.8(7 



TABLE I: Results for \Vub\, Jb, BR and the pull between BR 
and BR{B — ^ T^)exp in different scenarios (see text). 



using the UT fit and the result can be read from Eq. ([2| , 
Si?no-fit = (0.98 ±0.24) X 10-*. 

To summarize our results, we collect in Table |l] our 
predictions for BR in the considered scenarios. 



CONSTRAINTS ON THE MSSM PARAMETERS 



Bi?(B ^ ^i^)uuT = (3.9 ± 0.9) X IQ-^ (10) 
Si?(B ei^)uuT = (9-2 ± 2.1) x lO^^^ ^ 

where BR for these decays is defined in analogy with the 
B ^ TV case. 

Beyond MFV, the UUT construction is no longer ad- 
equate. Indeed, in the most general case, assuming only 
that NP contributions to semileptonic decays are negligi- 
ble, the prediction of BR{B — > tv) cannot be improved 



^ In the extraction of sin 2/3 from B — > Jj^K decays, possible 
NP enhancements of the penguin amplitude are bound using 
additional data l30l . 



It has been pointed out that the MSSM with MFV, 
TeV sparticles and large tan /3 could give negligible con- 
tributions to flavour physics except for B — )■ ri/, Amg, 
Bs A^^A^^ and B Xs^ 32J. We show that, with 
present data, the combination of the first three con- 
straints leaves little space for large tan /?. This can be 
easily understood as this model typically predicts a sup- 
pression of BR(B —J' Tv) rather than the enhancement 
required by the present measurements. An enhancement 
can be obtained only for very large values of tan /3 which, 
however, are disfavoured by the other constraints. 

We reanalyze the model of Ref. [32, with the follow- 
ing a-priori fiat ranges for the relevant low-energy SUSY 
parameters: \Ji = [-950, -450] U [450, 950] GeV, = 
[-3,3] TeV, tan/3 = [5,65], m^j+ = [100,1000] GeV, 
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FIG. 7: Regions in the (m^+,tan/?) parameter space of the 2HDM-II excluded at 95% probability by BR{B — >■ tv), BR{B — > 
Dtv)/BR{B -f Div) and BR{B -> X^-y). 



m, 



[400, 1000] GeV, m~g 



[400, 1000] GeV. The ex- 
pressions of i? — !■ TV, Bs — > ^-nd Arris can be found 
in Eqs. (3), (11) and (14) of Ref. |^ respectively. The ex- 
perimental constraints are Am^ = 17.77 ±0.12 ps~^ [33] 
and the upper bound BR{Bs — > < 5.8 x 10~® at 
95% C.L. [31]. 

In Figs. [s] we show the p.d.f. in the plane (tan/3, m^+) 
for /J, > 0. For completeness, in Figs. [9] and we present 
the corresponding one-dimensional p.d.f. for mfj+ and 
tan /?. As expected, the constraint from B ^ tv re- 
sembles the one obtained in the 2HDM analysis above 
(see Fig. [t]). Once the other constraints are included, 
however, the region at large tan j3/m[j+ is disfavoured. 
The combined exclusion region is roughly bounded by a 
straight line, giving tan/3 < 7.3m^+/(100 GeV) at 95% 
probability, with a remarkable similarity to the 2IIDM-11 
case. 

For < 0, the constraint from B ^ tv is less stringent 
for large tan/3, see Figs. [TT][r3] In fact, for /i < and 
very large tan /?, the interference with the SM in B 
TV becomes positive. However the combined bound is 



more severe than for /i > 0: for TO^f+ < 1 TeV, there 
is an absolute bound on tan/3 < 38 with at least 95% 
probability, while from the one-dimensional distribution 



in Fig. 13 we obtain tan/3 < 32 at 95% probability. 

For both signs of /j,, large values of tan /3 for sub- TeV 
charged Higgses are strongly disfavoured, including the 
fine-tuned region where the SUSY contribution enhances 
BR(B — >■ Tv) improving the agreement with the experi- 
mental average. 

From our analysis we also derive the following ranges 
for BR{Bs A^+M"): 



for ^ > 0, and 



[3, 8] X 10"^ @68% prob. 
[2, 26] X 10~^ @95% prob. 

[3, 6] X 10"^ @68% prob. 
[2, 17] X 10"^ @95% prob. 



(11) 



(12) 



for ^ < 0. These ranges can be compared with the SM 
prediction BR{Bs M+Ai")sM = (3.7 ±0.5) x 10"^. 
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FIG. 8: 68% (dark) and 95% (light) probability regions in 
the {mjj+ , tan j3) plane obtained using BR{B — )• tv) (top 
left), BR{Bs M^M") (top right), Am^ (bottom left), all 
constraints (bottom right) for /x > in the considered MFV- 
MSSM for the parameter ranges specified in the text. 



FIG. 10: 68% (dark) and 95% (hght) probability regions for 
tan/3 obtained using BR{B — >• tv) (top left), BR{Bs — >• 
fjb^ yL~) (top right), Arris (bottom left), all constraints (bot- 
tom right) for ^ > in the considered MFV-MSSM for the 
parameter ranges specified in the text. 
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FIG. 9: 68% (dark) and 95% (light) probability regions for 
mu+ obtained using BR{B — > ti/) (top left), BR{Bs — >■ 
IJb^ jjT) (top right), Ams (bottom left), all constraints (bot- 
tom right) for ^ > in the considered MFV-MSSM for the 
parameter ranges specified in the text. 



CONCLUSIONS 

We have shown how the use of the UT fit allows to im- 
prove the prediction of BR{B tv) in the SM, thanks 
to a better determination of \Vub\ and /s. Considering 
the generalization of the UT fit to various NP scenarios, 



FIG. 11: Same as Fig.[8]for /i < 0. 

we have obtained results for BR^ defined as the predic- 
tion of BR{B — )• Tv) assuming negligible NP contribu- 
tions to the decay amplitude. The comparison of BR 
to the experimental result provides an improved probe 
of the presence of NP in the decay amplitude. Our re- 
sults are summarized in Table [ij Finally, we studied the 
present constraints on the 2HDM-II and on the MFV- 
MSSM with TeV sparticles. In both models, we find that 
large values of tan f3 for sub-TeV charged Higgs masses 
are disfavoured by present data. 
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FIG. 12: Same as Fig.[9]for fi< 0. 




tanp tanp 

FIG. 13: Same as Fig.[lO]for < 0. 

We thank G. Isidori for useful comments and discus- 
sion on our analysis. We acknowledge partial support 
from RTN European contracts MRTN-CT-2006-035482 
"FLAVIAnet" and MRTN-CT-2006-035505 "Heptools". 
M.C. is associated to the Dipartimento di Fisica, Univer- 
sita di Roma Tre. E.F. and L.S. are associated to the Di- 
partimento di Fisica, Universita di Roma "La Sapienza" . 



[1] S. L. Glashow, J. Iliopoulos and L. Maiani, Phys. Rev. D 
2 (1970) 1285. 

[2] B. Grzadkowski and W. S. Hon, Phys. Lett. B 272 (1991) 



383; Phys. Lett. B 283 (1992) 427. 
[3] C. Amsler et al. [Particle Data Group], Phys. Lett. B 667 
(2008) 1. 

[4] V. Lubicz and C. Tarantino, Nuovo Cim. 123B (2008) 

674 arXiv:0807.4605 [hep-lat]]. 
[5] N. Cabibbo, Phys. Rev. Lett. 10 (1963) 531; 

M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49 

(1973) 652. 

[6] T. Hokuue et al. [Belle Col laboration] , Phys. Lett. 

B 648 (2007) 139 [arXiv:hep -ex/0604024 ; B. Aubert 

et al. [BABAR Collaboration], Phys. Rev. Lett. 97 

(2006) 211801 arXiv:hep-ex/0607089 ; Phys. Rev. Lett. 

98 (2007) 091801 arXiv:hep-ex/0612020 ; Phys. Rev. 

Lett. 101 (2008) 081801 arXiv:0805.2408 [hep-ex]]; 

N. E. Adam et al. [CLEO CoUaboration] , Phys. Rev. 

Lett. 99 (2007) 041802 arXiv:hep-ex/0703041 ; I. Adachi 

et al. [Belle Collaboration], arXiv:0812.1414 [hep-ex]. 
[7] E. Barberio et al. [Heavy Flavor Averaging Group], 

arXiv:0808.1297 [hep-ex]. 
[8] P. BaU and R. Zwicky, Phys. Rev. D 71 (2005) 014015 

I arXiv:hep-ph/04 06232. . 
[9] B. O. Lange, M. Neulaert and G. Paz , Phys. Rev. D 

72 (2005) 073006 arXiv:hep-ph/0504071 ; J. R. Ander- 

sen and E. Gardi, JHEP 0601 (2006) 097 |arXiv:hep- 

I ph/0509360; U. Aghetti, G. Ferrera and G. Ric- 

dardi, Nucl. Phys. B 768 (2007) 85 |arXiv:hep- 

ph/0608047; P. Gambino et al, JHEP 0710 (2007) 058 

[arXiv:0707.2493 [hep-ph]]; U. Aghetti et al., Eur. Phys. 

J. C 59 (2009) 831 arXiv:071 1.0860 [hep-ph]]; E. Gardi, 

'arXiv:0806.4524 [hep-ph]. 
[10] K. Ikado et al. [Belle Collaboration], Phys. Rev. Lett. 97 

(2006) 251802 arXiv:hep-ex/0604018 ; B. Aubert et al. 
BABAR Collaboration], Phys. Rev. D 77 (2008) 011107 
arXiv: 0708.2260 [hep-ex]]; arXiv:0809.4027 [hep-exh 
. Adachi et al. [Belle Collaboration], arXiv:0809.3834l 

[hep-ex] . 

[11] M. Bona et al, arXiv:0709.0451 [hep-ex]. 

[12] M. Bon a et al. [UTfit Collaboration], JHEP 0610 (2006) 

081 [arXiv:hep-p h/0606167 . 
[13] M. Bona et a l. [iJTfit Collaboration], JHEP 0507 (2005) 

028 [arXiv:h ep-ph/0501199. 
[14] J. Charles et al. [CKMfitter Group], Eur. Phys. J. C 41 

(2005) 1 arXiv:hep-ph/0406184 . 
[15] M. Artuso et al. [Cleo Collaboration], Phys. Rev. Lett. 

75 (1995) 785; N. Satoyama et al. [Belle Collaboration] , 

Phys. Lett. B 647 (2007) 67 |arXiv:hep-ex/0611045i ; 

B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 

77 (2008) 091104 arXiv:0801.0697 [hep-ex]]. 
[16] B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 

79 (2009) 091101 arXiv:0903.1220 [hep-ex]]. 
[17] R. S. Chivukula and H. Georgi, Phys. Lett. B 188 (1987) 

99; L. J. HaU and L. RandaU, Phys. Rev. Lett. 65 (1990) 

2939; E. Gabrielh and G. F. Giudice, Nucl. Phys. B 433 

(1995) 3 [Erratum-ibid. B 507 (1997) 549] |arXiv:hep- 

I lat/9407029; M. Misiak, S. Pokorski and J. Rosiek, 

Adv. Ser. Direct. High Energy Phys. 15 (1998) 795 

|arXiv:hep-ph/9703442 ; M. Ciuchini et al., Nucl. Phys. B 

534 (1998) 3 arXiv:hep-ph/9806308 ; A. J. Buras, Acta 

Phys. Polon. B 34 (2003) 5615 arXiv:hep-ph/0310208;;; 

M. Blanke et al., JHEP 0610 (2006) 003 |arXiv:hep^ 
' ph/0604057. 

[18] G. DAmbrosio et al, Nucl. Phys. B 645 (2002) 155 

|arXiv:hep-ph/0207036 ; 
[19] M. Bona et al. [UTfit Collaboration], JHEP 0803 (2008) 



9 



049 arXiv:0707.0636 [hep-ph]]. 
[20] M. Bona et al. [UTfit Collaboration], |arXiv:0803.0659 
[hep-ph]. 

[21] A. J. Buras et al, P hys. Lett. B 500 (2001) 161 
(arXiv:hep-ph/0007085| . 

[22] W. S. Hou; Phys. Rev. D 48 (1993) 2342. 

[23] A. Heister et al. [ALEPH Collaboration], Phys. Lett. 
B 543 (2002) 1 arXiv:hep-ex/0207054 ; J. Abdallah et 
al. [DELPHI Collaboration], Eur. Phys. J. C 34 (2004) 
399 arXiv:hep-ex/0404012 ; P. Achard et al. [L3 Col- 

laboration], Phys. Lett. B 575 (2003) 208 |arXiv:hep-| 

ex/0309056 . 

[24] M. Misiak et al, Phys. Rev. Lett. 98 (2007) 022002 

[arXi v:hep-ph/0 609232 . 
[25] M. T anaka, Z. Phys. C 67 (1995) 321 arXiv:hep- 
ph/9411405 ; H. Itoh, S. Komine and Y. Okada, Prog. 

Theor. Phys. 114 (2005) 179 arXiv:hep-ph/0409228 . 
[26] B. Aubert et al. [BABAR CoUaboration] , Phys. Rev. 

Lett. 100 (2008) 021801 arXiv:0709.1698 [hep-ex]]; 



Phys. Rev. D 79, 092002 (2009) [arXiv:0902.2660 [hep- 
ex]]; T. lijima, talk given at the 24th International Sym- 
posium On Lepton-Photon Interactions At High Energy 
(LP09), http : //tinyurl . com/ii j imaLPOQ. 

[27] J. F. Kamenik and F. Mescia, Phys. Rev. D 78 (2008) 
014003 |arXiv:0802.3790 [hcp-ph]].' 

[28] H. E. Haber, arXiv:hep-ph/9707213" 

[29] O. Deschamps et al, arXiv:0907.5135 [hep-ph]. 

[30] M. Ciuchini, M. Pierini and L. Silvestrini, Phys. Rev. 
Lett. 95 (2005) 221804 |arXiv:hep-ph/050729(|^ 

[31] M. Ciuchini et al., JHEP 0107 (2001) 013 arXivihep- 
r ph/0012308 . 

[32] G. Isidori and P. Paradisi, Phys. Lett. B 639 (2006) 499 

|arXiv:hep-ph/0605012 . 
[33] A. Abulencia et al. [CDF Collaboration], Phys. Rev. Lett. 

97 (2006) 242003 arXiv:hep-ex/0609040 . 
[34] T. Aaltonen et al. [CDF CoUaboration], Phys. Rev. Lett. 

100 (2008) 101802 |arXiv:0712.1708| [hep-ex]]. 



